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Abstract. Formation and disruption of dark-matter cusps are reviewed. Accumulation 
of baryons at the center of a halo can displace the dark matter, converting singular 
density cusps into low-density cores. The displaced mass can be of order ~ 10M. with 
M, the mass of the infalling population. If M. is identified with the masses of the 
black holes currently observed at the centers of bright galaxies, predicted core radii 
are ~ a few x 10 2 pc. Other mechanisms, such as early mass outflow, may explain the 
large dark-matter cores in dwarf and low-surface-brightness galaxies. Predictions of 
dark matter annihilation radiation from the center of the Milky Way galaxy are shown 
to be strongly dependent on the galaxy's merger history. 

1 CDM Cusps 

The influential paper of Navarro, Frenk & White jlj proposed a universal, broken 
power-law profile for dark matter halos in cold dark matter (CDM) cosmologies. 
The dependence of p on distance from the halo center was found to be well 
described by a two-parameter family of distributions, 

PNF\v(r) = PoC 1 (1 +C)~ 2 , i = r/r a . 

The central density of the NFW profile diverges as p ~ r , similar to what was 
found in earlier iV-body studies (e.g. Q). However the small-radius dependence 
was little more than an ansatz since the relevant scales were barely resolved in 
the simulations. A debate ensued as to whether the profiles arc indeed universal, 
and if so, what power of radius describes the dark matter density in the limit 
r — > 0. A number of subsequent studies (^,§|||||j7j found steeper central profiles, 
e.g. p ~ r — 1-5 , the "Moore" profile. 

At about the same time, the luminous (stellar) densities in galaxies were 
found to be centrally divergent ||||; in contrast with CDM halos, the power law 
index of stellar nuclei varies widely but correlates well with galaxy properties, 
in the sense that fainter galaxies have steeper cusps JlO| . While the power-law 
nature of stellar cusps suggests a gravitational origin, gas dynamics probably also 
played a role in their formation, including processes like radiative dissipation, 
star formation, and energy infusion from supernova explosions. 

The power-law nature of density stratification in virialized objects arises 
naturally in the standard cosmological paradigm where the growth of objects 
proceeds from smallest to largest scales starting from a featureless spec- 
trum of primordial fluctuations P(k) oc k n . Building blocks of structure can 
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be identified with peaks in the field of primordial Gaussian fluctuations. The 
ensemble-averaged "profile" of an individual density peak is then a power law 
with slope related to the tilt n of the primordial spectrum on a given scale 
[ ^2p^ , p^ |. Furthermore, evolution of virialized systems through mergers tends 
to preserve the power-law slope of the central cusp |H| . Outside of the central 
cusp, the profile will be modified by the tidal field of neighboring peaks and by 
the redistribution of material in unequal-mass mergers. 

It was noticed rather quickly following the advent of realistic large-scale cos- 
mological A^-body simulations that after the slope and normalization of the den- 
sity profiles are taken into account, cusps predicted by the simulations are too 
concentrated to accommodate dynamical constraints derived from the rotation 
curves of (especially low surface brightness) galaxies jL(|[l7|. Better fits were 
obtained with ad hoc, constant-central-density profiles, e.g. the non-singular 
"isothermal sphere" Core sizes are not well constrained but recent 

model fits to the rotation curve data give r c < 1 kpc for most galaxies (e.g. 
po|). Other observations on galaxy scales (such as strong lensing |^l|,^2)) have 
led to more equivocal assessments of the compatibility of dark matter models 
with observations. Meanwhile, the most detailed study to date of resolution ef- 
fects in the CDM simulations |23] argues that the innermost resolved cusp slope 
is marginally steeper than NFW, p ~ r -12 , and that it exhibits a shallowing 
trend with decreasing radius. 

The cusp problem in CDM cosmologies can be alleviated either by postulating 
changes in the basic physics, or by adding baryons that couple to the dark matter 
and modify its profile. Proposed changes in fundamental physics include weakly- 
interacting particles with large scattering cross-sections | j24p5|| , "warm" dark 
matter |2(|, broken scale-invariance (27J, or modified gravity [ [28|]29| ]. Standard- 
physics solutions include outflow of gas during an early phase of galaxy formation 
[ p0| , torques from a barlike perturbation J3l]j2| or infall of a population of 
condensed objects which disrupt the cusp fl33| , M ■ 

Current notions of galaxy formation are not precise enough to allow testing 
of most of these ideas. Here we focus on the last mechanism, displacement of 
dark matter by condensed objects, because the dynamics are straightforward and 
the predicted level of the effect is quantifiable. A lower limit on the displaced 
mass is set by the known masses of compact objects at the centers of galaxies, 
the supermassive black holes (SBHs) j$5), some of which are known to have 
been present since z ~ 6 Regardless of their detailed formation history, 
these black holes must have grown partly through the coalescence of less-massive 
black holes during galaxy mergers, displacing and ejecting dark (and baryonic) 
matter. Predicted core sizes are ~ hundreds of parsecs, consistent with the sizes 
of stellar cores in moderate- to high-luminosity galaxies, which probably formed 
in the same way. While it remains to be seen whether this model can explain 
all of the discrepancies between predicted and observed densities at the centers 
of halos, particularly in low-surface-brightness and dwarf galaxies, the model is 
almost certainly relevant to the dark matter distribution within a few tens of 
parsecs of the SBHs where dark matter annihilation radiation would be peaked. 
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2 Enhancement and Disruption of Cusps 

Accumulation of baryons at the center of a dark-matter halo will deepen the grav- 
itational potential and steepen the cusp, but the net effect depends on whether 
the baryons also deposit energy into the dark matter as they fall in. The simplest 
model (Fig, la) is adiabatic growth of a central point mass due to radial infall 
of gas [p8|j37| . Let the initial mass of the central lump be zero and its final mass 
Mm- For simplicity, assume that the orbits of the dark matter particles are cir- 
cles; if growth of the central lump is slow, the orbits will remain circular due to 
adiabatic invariance. Conservation of angular momentum implies riMi — rfM, 
where fj and rt are the initial and final radii of a dark-matter particle, Mi is 
the dark mass within and M, is assumed to dominate the potential within r f . 
Conservation of dark matter further implies Mi{r,i) = M/(r/). 

If the initial dark matter density profile is a power law, pi ~ r -7 , the two 
conservation laws imply pj ~ r~ s , S — (9 — 2j)/(A — 7). A more realistic, 
isotropic orbital distribution predicts the same relation between 7 and S f^fl . 
For < 7 < 2, the index of the power-law "spike" so formed lies between 2.25 
and 2.5; the spike extends roughly to the outer radius of influence of the central 
object, r ~ rft = GM./er 2 , where a is the velocity dispersion of the dark-matter 
particles in the pre-existing halo. In physical units, 

r h w 10M 8 ct^oo pc 

where the central mass has been normalized to 10 8 M Q and the velocity dispersion 
to 200 km s _1 . This model has been widely used to predict the flux of annihilation 
radiation from dark matter particles around the Milky Way black hole; most of 
the radiation would come from the spike [ ^9]j40|] . 

The outcome is very different if the lumps form off-center, or are deposited in 
the halo following a merger (Fig. lb). The lumps spiral inward due to dynamical 
friction, transferring their energy into the background and lowering the central 
density H||P,||,||. One well-studied case is the formation of a binary SBH 
following a merger. The binary first forms at a separation of a ~ GM m /a 2 , 
with Mm = mi + 7Ti2 the binary mass. Continued evolution of the binary occurs 
through the gravitational slingshot : dark matter particles that approach 
the binary within a distance ~ 3a are kicked out at velocities of ~ (GM./a) 1 / 2 . 

Slingshot ejection occurs at a rate determined by the coupled equations 

dt\a) a ' 

dM e ~ , . 

dUYfa) - JM ' (2) 

| ft5| . Here M e j is the mass ejected by the binary and p and a refer to the back- 
ground. The dimensionless numbers H and J are functions of mi/m-2 and the 
binary separartion; in the limit of a "hard" binary, a -C rj ly and for mi ~ m 2 , 



H ~ 16 and J ~ 1 46 1. Solving the coupled equations for M e j(t) is not straight 



forward however since p changes with time as the dark matter is ejected. Simple 
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Fig. 1. Two schemes for accumulating baryons at the center of a dark- matter halo, (a) 
Mass falls in radially without depositing energy into the dark matter. The central po- 
tential deepens and the cusp responds by becoming steeper, (b) Lumps form off-center 
and spiral in, as they would following a merger. As their orbits decay via dynamical 
friction, the lumps transfer energy to the dark matter and lower its density. 



models that take into account this time dependence (47]] give M e j cx M, In t, 
consistent with what is seen in detailed iV-body simulations [4q|. 



Fig. 2 shows the results of one such simulation 1 48 . The merger was between 
equal-mass spheroids with steep central cusps, p ~ r~ 2 , and point particles 
representing the two SBHs. The r~ 2 cusp is gradually converted into an r _1 
cusp by the gravitational slingshot; at the end of the simulation, the ejected 
mass is roughly twice the mass of the binary SBH. 

The total mass ejected by a point-mass binary depends on how long the 
binary continues to shrink. An upper limit in the case of binary SBHs is given 
by the mass ejected in decaying to a separation of a gr , the distance at which 
emission of gravitational radiation induces a rapid coalescence. The ejected mass 
is 

M ej *m 2 ln(^-) (3) 



with ni2 the mass of the smaller black hole and ah ~ GM.jo 2 47 . For binary 
SBHs, Mm ~ 10 8 A/ Q , a gr ~ 0.01 pc and the total mass ejected in decaying to a gr 
is ~ 5M,. However a binary may "stall" after ejecting less mass if the density 
in its vicinity drops too low 

If the number of massive infalling objects exceeds two, heating of the back- 
ground by binaries can occur indirectly as well: massive particles are given kicks, 
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Fig. 2. Formation of a core by a binary SBH |4q ]. (a) Density profiles. Dashed profile 
is the initial model, p ~ r -2 ; solid lines show the evolution of p starting roughly at the 
time when a hard binary is formed, (b) Surface density profiles. 



then lose their kinetic energy to the background as they spiral in and are kicked 
out again etc. As a crude model of this process, consider the evolution of an 
ensemble of point particles of mass m, orbiting in a dark-matter halo with 
p ~ r~ 2 - slightly steeper than the cusps found in CDM simulations but eas- 
ier to treat. The energy released as one lump spirals in from radius to is 
2m. a 1 ln(rj/r/). Decay will halt when the lumps form a self-gravitating system 
of radius ~ GM./c 2 with M, the total baryonic mass. Equating the energy re- 
leased during infall with the energy of dark matter initially within r c , the "core 
radius," gives 

2GM, fr t a 2 \ 

Most of the particles that deposit their energy within r c will come from radii 
j-j ss a few x r c , implying r c ps several x GM, / a 2 and a displaced mass of 
~ several x M,. 

Evolution will continue as the massive particles form binaries and begin to 
kick out other massive particles. Assume that ejection occurs via many small 
kicks, such that almost all of the binding energy so released can find its way into 
the dark matter as the ejected particles spiral back into the core. The energy 
released by a single binary in shrinking to a separation such that its orbital 
velocity equals the escape velocity from the core is ~ m.cr 2 \n(AM,/Mhaio)- If 
all of the massive particles find themselves in such binaries before their final 
ejection and if most of their energy is deposited near the center of the halo, the 
core mass will be 

M e ,«M.ln(^) (5) 
e. g. - 10M. for M./M halo « lfr 4 . 
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Fig. 3. Formation of a core by infalling lumps; details of the simulation are given in 
the text, (a) Initial and final density profiles of the dark matter, (b) Ejected dark mass 
as a function of time, defined as the difference in enclosed mass between the profiles 
at times and T inward of their first intersection. M, is the total "baryon" mass, i.e. 
the mass of the infalling population. T c is the crossing time. Almost all of the massive 
objects have been ejected by the end of the simulation. 

A detailed A^-body simulation illustrates this process (Fig. 3). The initial 
model for the dark-matter halo has p ~ r -15 at small radii and total mass M. 
The condensed objects were represented by 20 point masses, each weighing 10 
times as much as a dark matter particle and with combined mass M, = 10~ 2 M. 
Evolution occurs in two stages, as discussed above. As the orbits of the massive 
particles initially decay, they transfer energy to the background and lower the 
central density. Once the massive particles form a bound subsystem, further 
evolution is driven by the formation and disruption of binaries. The result is 
a slower growth of the dark matter core, as massive particles arc ejected from 
the central regions and fall back in. At the end of the simulation the dark mass 
ejected from the cusp is ~ SM,. 

Only a single massive binary is left at the end - most of the particles that 
displaced the dark matter have been ejected from the system. This model is the 
dynamical analog of wind-driven outflow models |30|pH : in both schemes, the 
baryons that created the core are no longer in evidence after the core is in place. 

Is this a realistic model? Multiple baryonic clumps are often postulated to 
form via gaseous dissipation in the early stages of galaxy formation (e.g. ]3l| , |3^] ) 
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Fig. 4. Cusp destruction in hierarchical mergers of halos containing central "black 
holes" El|. Left panel: a merger tree; time proceeds downwards and the area of each 
circle is proportional to the mass of the halo, and of its central SBH. Right panel: Halo 
profiles. The top profile is the initial halo model (e.g. U, V); the middle profile is halo 
T; the bottom profile is halo S. The profiles are scaled so that the black hole mass is 
the same in each. The progenitor halos were truncated Moore profiles, with p ~ r -1,5 
central density cusps. Axis scalings are based on an assumed black hole mass of 1O 8 M , 
a halo mass of 1O 13 M0 and a virial radius of 300 kpc. 



but their masses, numbers and degree of central concentration are poorly known. 
Another possible candidate for the lumps is an early population of black holes, 
e.g. the remnants of Population III stars (to. ~ 150M Q ) pofl . But current theory 
suggests that only one such black hole forms per mini-halo at a redshift of ~ 20 
plj , and furthermore the inspiral time for an object of this mass would be of 
order 10 9 yr, unless it came attached to a larger mass. Black holes of intermediate 
mass (to. ~ 1O 3 M0) are a third possibility but their existence is debated and 
their time of formation uncertain. 

SBHs have been present in at least some spheroids since a redshift of ~ 6 
[p6| ; many of the SBHs that we see now should therefore be the products of 
multiple mergers. Typically those mergers will have been "minor," to-2 <C mi. 
In such an interaction the ejected mass scales with TO2, the mass of the smaller 
hole; the total mass ejected in a sequence of such mergers would then depend 
only on Af., the final mass. However if a black hole grows via a merger hierarchy 
of comparably-massive black holes, the effect on the density cusp should be to a 
certain extent cumulative: a merger of two galaxies whose cusps had previously 
been destroyed by binary black holes, will produce a shallower central profile than 
a merger between two galaxies with initially steep cusps. This is verified via an 
A^-body experiment in Fig. 4. The ejected mass scales both with M., the final 
black hole mass, and N, the number of stages in the merger hierarchy. 
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There is some quantitative empirical support for idea that cores are produced 
from the binding energy released by binary black holes. The luminous (stellar) 
matter in elliptical galaxies and bulges exhibits steep power-law cusps in faint 
galaxies and cores in bright galaxies among the "core" galaxies, the core 
mass is of order ~ 10M. with M, the current black hole mass [Q. The largest 
cores observed, in galaxies like M87, have sizes of ~ 0.5 kpc; it seems circum- 
stantially reasonable to suppose that dark-matter cores could be comparable in 
size. The persistence of steep cusps in faint elliptical galaxies is puzzling but 
may indicate that the stellar cusps formed by dissipative processes during or 
after the last significant merger |53|] ; in any case, presence of a stellar cusp is not 
necessarily evidence for a dark matter cusp. 

SBHs represent an extreme in terms of compactness, but galaxies contain 
other condensed components that could in principle act like "black holes" during 
mergers, displacing dark matter as they spiral inward. Bulges are an obvious 
candidate; with masses ~ 10 3 times those of SBHs Q they could in principle 
displace a much larger quantity of dark matter. In order for two bulges to act like 
a "binary black hole," they must be compact enough to not overlap appreciably 
at the separation where they first form a bound subsystem. This separation 
is ~ GMbuige/a 2 with a the dark matter velocity dispersion. Requiring this 
separation to be greater than the bulge scale length Rbuige gives GMbuige 3> 
RbuigeV 2 or Ubuige cr, with Obuige the bulge velocity dispersion. Thus a "hot" 
bulge is a "hard" bulge. In fact Obuige appears to be comparable to, though 
somewhat larger than, a J5j| so it is unlikely that bulges - at least as we currently 
observe them - would be effective at displacing dark matter during mergers. 

3 An Application: Annihilation Radiation from the 
Galactic Center 

If dark matter consists of massive weakly interacting particles like neutralinos, 
the particles should self-annihilate producing a potentially observable signal. The 
annihilation rate in a spherical halo is 

F= W* /ArH^dr, (6) 

with M x the ncutralino mass, av the annihilation cross section, and p the dark 
matter density [^6| . Annihilation products include gamma rays, from the decay 
of neutral pions; high-energy neutrinos; and synchrotron emission from high 
energy positrons and electrons moving in the ambient magnetic field. Since the 
annihilation flux depends on the squared density of neutralinos, the signal would 
be greatly enhanced along lines of sight where the dark matter is clumped. These 
include the center of the Milky Way galaxy, where the density in a smooth halo 
would be maximum, as well as directions that intersect the centers of relic halos 
orbiting as subclumps in the Milky Way halo |57| , p8| . 

The annihilation signal from the galactic center would be further enhanced if 
there is a power-law spike associated with the central SBH, as discussed above. 
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Fig. 5. Dark matter annihilation radiation from the galactic center |62|]. Left panel: 
Density profile of the Milky Way dark halo. Curves labelled M and S show the halo 
profile before and after growth of the central SBH. Other curves show profiles following 
mergers with smaller halos with various mass ratios: (A, B) 1:1; (C, D) 1:3; E 1:5; F 
1:10. The mass of the SBH in the smaller halo scales with halo mass. Right panel: 
Dimensionless integrated annihilation flux as a function of the angular acceptance of 
the detector, centered on the SBH. 



Gondolo & Silk |39[ considered this case and showed that the flux is increased 
by ~ 5 orders of magnitude compared with models lacking a central spike. This 
prediction has been claimed to be inconsistent with measured fluxes of annihi- 
lation products, including synchrotron emission Q and neutrinos Q (but see 

S3). 

The central slope of the Milky Way's dark matter distribution could have 
been substantially lowered if the halo experienced a merger after formation of 
the spike, assuming that the infalling halo also contained a SBH. At least one 
such merger, with a mass ratio of 1 : 10 or less, is statistically likely to have been 
experienced by a halo as massive as that of the Milky Way since a redshift of 
2 |62| . Fig. 4 shows density profiles of merged dark matter halos, in simulations 
where the initial profile of the larger halo was normalized to match the Milky 
Way, with a central r~ 15 density cusp J6^]. The spike is efficiently destroyed by 
the binary SBH in all of the mergers, greatly reduced the predicted annihilation 
flux. 

Cusp destruction may have been avoided in the case of our galaxy ]6l|| . Little 
is known about the detailed merger history of the Milky Way but it has been 
argued that the ages of stars in the galaxy's "thick disk" imply a time of ~ 11 Gyr 
since the last significant merger |63|] . The Milky Way SBH may have acquired 
most of its mass during or even after this merger event. However, if the SBH 
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grew significantly after a merger that had destroyed a pre-existing dark matter 
spike, it is unlikely that the new spike would be as steep. The reason is that 
the dark matter density would probably be non-singular following the merger 
(e.g. Fig. 4), and growth of a point mass in a non-singular background produces 
a shallow spike, p ~ r -1 ' 5 [ |37| . The annihilation flux from such a spike would 
be only slightly greater than the expected background in the direction of the 
galactic center p9[ . 
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